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Abstract 

A breadth-first search method for determining optimal 3-line circuits composed of quantum NOT, CNOT, 
controlled-V and controlled-V+ (NCV) gates is introduced. Results are presented for simple gate count and 
for technology motivated cost metrics. The optimal NCV circuits are also compared to NCV circuits derived 
from optimal NOT, CNOT and Toffoli (NCT) gate circuits. The work presented here provides basic results and 
motivation for continued study of the direct synthesis of NCV circuits, and establishes relations between function 
realizations in different circuit cost metrics. 

1 Introduction 

Reversible and quantum logic synthesis have attracted recent attention as a result of advances in quantum and nano 
technologies. Many of the proposed synthesis methods, especially in the area of reversible logic synthesis, assume 
large gate libraries where implementation costs of the gates may vary significantly. However, these methods most 
often target minimization of the gate count. Use of such a circuit cost metric is Ukely to result in seemingly small 
circuits, which are in fact expensive to construct. 

The synthesis of circuits composed of NOT, CNOT and Toffoli (NCT) gates 1141 and multiple control Toffoli 
gates ^9, jl6, ^ li, i4J has recently been extensively studied. Converting an NCT or a multiple control Toffoli gate 
circuit into one composed of NOT, CNOT, controlled-V and controlled-y+ (NCV) gates has also been considered 
||2]E1. Further simplification of NCV circuits has also been well studied Q- 

It is an important question how close the circuits found by the above approaches are to optimal. The direct 
synthesis of NCV circuits is also of considerable interest since intuitively one would expect that to produce better 
results than the indirect route via NCT circuits. Finally, we note that observing optimal circuits for small cases 
often will shed light on good (if not optimal) synthesis approaches. 

For these reasons, we here present an approach to finding optimal 3-line NCV circuits using various cost 
metrics. We compare these results to those found through mapping optimal NCT circuits. The advantage of direct 
NCV synthesis will be clear even for the 3-line case. Also, the results clearly demonstrate the difference between 
using simple gate count and technology motivated cost metrics. 

The necessary background is reviewed in Section |2] A breadth-first search procedure to find optimal 3-line 
NCV circuits is given in Section |3l and properties identified in those circuits are discussed in Section |3 Section 
|5] presents comparative results to circuits derived from NCT circuits and for various cost metrics. The paper 
concludes with remarks and suggestions for ongoing research in Section^ 



2 Background 

We here provide a brief review of the basic concepts required for this paper. For a more detailed and formal 
introduction we refer the reader to [11]. 
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A single quantum bit (qubit) has two values, or 1, traditionally depicted as |0) and 1 1) respectively. The state 
of a single qubit is a linear combination a|0) + P|l) (also written as a vector (a, P)) in the basis {|0), 1 1)}, where 
a and p are complex numbers called the amplitudes, and |ap + |pp = 1. Real numbers |ap and |pp represent 
the probabihties (p and q) of reading the values |0) and |1) upon physical measurement of the qubit. The state 
of a quantum system with « > 1 qubits is described as an element of the tensor product of the single state spaces 
yielding a normalized vector of length 2" called the state vector Quantum system evolution results in changes of 
the state vector expressible as products of 2" x 2" unitary matrices. This formulation characterizes a transformation 
but provides no indication of its implementation cost. 

In one common approach, small gates are used as elementary building blocks with unit cost 1111 121131171151. 
Among them are NOT {x x) and CNOT {{x,y) — *■ {x,x(By)) gates, the 2-bit controlled-V gate which changes 

the target line according to the transformation given by matrix V = ^ ^ ^ if, and only if, the control line is 1 . 

Controlled- V+ analogously applies the transformation V+ = V^'. Consult pi H for the case where the control line 
assumes a quantum (non-Boolean) value. Gates controlled-V and controlled-y+ are also referred to as controlled- 
sqrt-of-NOT gates since — ^NOT. In order to implement a Boolean specification, and assuming an 

auxiliary line is available in addition to the minimal set of lines needed for reversibility |6|, the set of gates NOT, 
CNOT, controlled-y and controlled-y+ is complete. We call this set the NCV gates. 

Definition 1. The NCV-111 cost of a circuit composed of NCV gates is the number of gates in the circuit. 

In an alternative approach II0III3I . it is observed that any circuit can be composed with single qubit and CNOT 
gates and the circuit cost is calculated based on the number of CNOT gates required. With regards to the NCV 
library, we next define NCV-012 cost and motivate our definition by the fact that controlled-V and controlled-y+ 
gates require at most 2 CNOT gates when decomposed into a circuit with single qubit and CNOT gates ( II II . page 
181). 

Definition 2. The NCV-012 cost of an NCV circuit is linear with weights 0, 1, 2, and 2 associated with the gates 
NOT CNOT controlled-V and controlled-V^ , respectively. 

The last cost metric that we consider in this paper is motivated by a recent investigation of the technological 
(liquid NMR) costs of quantum and reversible primitives Js] . 

Definition 3. The NCV-155 cost of an NCV circuit is linear with weights 1, 5, 5, and 5 associated with the gates 
NOT CNOT controlled-V and controlled-V^ , respectively. 

Finally, we define the multiple control Toffoli gate |15l. 

Definition 4. For the set of Boolean variables {x\,X2,...^ x„} the Toffoli gate has the form TOF{C;T), where 
C = {x,| ,x/2, ...,JC,j}, T = {xj} and CCiT = 0. It maps the Boolean pattern {x'^ ,X2 , ■■■,x'^) to (x^^jXj", , 
xj ®x^xJ...jCj^,xt^j, ...,x^). C will be called the control set and T will be called the target. 

The Toffoli gate and its generalizations with more than two controls form a good basis for synthesis purposes 
and have been used by many authors ||9][l6||8|[l]@l. However, Toffoli gates are not simple transformations. Rather 
they are composite gates themselves and Toffoli gates with a large set of controls can be quite expensive ||2|E1. As 
a result, the NCV- 11 1 cost of a 100-gate Toffoli circuit with 10 input/output lines can be as low as 100, or may be 
as high as 15,200. These are, of course, extreme numbers, however, in Section|5]we show a better analysis of how 
the costs may differ. 

3 The Search Procedure 

Our main tool for the investigation of relations amongst optimal NCV-111, NCV-012 and NCV-155 implemen- 
tations and their relation to the most commonly used reversible circuit cost metric, multiple control Toffoli gate 
count, is a search procedure for the optimal synthesis of NCV circuits for all 40,320 size 3 reversible Boolean 
functions. For this problem, we use a prioritized, pruned breadth-first search. This is because the search tree for 
simple full breadth-first search grows too fast for a search to be accomplished in feasible time and space. This is 
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in contrast to using a full breadth-first search for NCT optimal synthesis 1141 . which is tractable for 3-line circuits 
with no search optimizations. 

To see the size of the problem, note that the number of base transformations for NCV gates is 21: 3 trans- 
formations involve a NOT gate, and the use of the gates CNOT, controlled-V and controlled-V^ each result in 6 
transformations. The results in Table|2]show that the length of the optimal implementations can be up to 16. Thus, 
the number of nodes at one level of the search tree can be as high as 21'^ « 1.4* 10^'. 

Breadth first search was previously applied to the synthesis of optimal circuits for the size 3 reversible functions 
usmg NOT, CNOT, and Toffoli gate library |I4|. The size of the bottom level of their search tree (branching factor 
to the power of tree depth) is 9*^ ~ 4.3 * 10^ showing that such a search is significantly simpler than the one we are 
pursuing and requires no techniques to reduce the search space or make the search efficient. In addition, Toffoli 
gates are not simple transformations, while synthesis of optimal circuits makes more sense in terms of simple 
transformations. Our search procedure can find optimal circuits in any weighted gate count metric, not just in the 
simple gate count metric used in 1141 . 

An earlier attempt to synthesize optimal NCV circuits |I17| was capable of synthesizing optimal implementa- 
tions of the maximal cost 7 in one specific cost metric (NCV-011?). The maximal number of synthesized functions 
is 10136 1171 . which is about a quarter of the pool size. Our program synthesizes optimal implementations for 
all 40,320 size 3 reversible functions, and is not tied to a specific costing metric. Further, our search procedure is 
significantly faster: it completes the entire search in the lesser time than the search of all optimal 5-gate implemen- 
tations in 1171 . 

In our work, we do not allow a quantum gate (controlled-V and controlled-y+) to have a control line that may 
at that point in the circuit take a quantum (non-Boolean) value. We do not have a mathematical proof that deleting 
this restriction will not result in construction of smaller circuits. However, numerous experiments indicate that use 
of a quantum gate with a quantum control line does not lead to a more efficient circuit than the one we construct. 
The above restriction allows us to work with quaternary logic instead of a continuous values. 

The techniques used to reduce the size of the search include: 

• Based on the observation above, we can view a quantum function as a base reversible Boolean function with 
a quantum factor added. In other words, each of the values that may occur in the truth table is stored as a 
2-bit number Values 0, l,V and y+ are represented by 00, 10, 01 and 11, respectively. 

• While our final goal is to synthesize all size three reversible functions, and none of the quantum (non- 
Boolean) functions, the latter still have to be stored and referred to during the search process. Each such 
function is stored in a queue associated with its base reversible Boolean function that can be uniquely identi- 
fied from the first digit in the above 2-bit encoding. The quantum function itself is then identified by a 24-bit 
number composed of the second bits, that we call the quantum signature. A quantum function is uniquely 
defined by the base Boolean function and its quantum signature. 

• When we assign a new gate to the existing optimal cascade we never choose a gate with the same set of 
controls and targets as the immediately previous one in the circuit. This is because such a sequence can 
always be reduced by template application |7j, and thus will not be part of an optimal circuit. This cuts 
down the number of gates that one must consider at each step. 

• We note that in an NCV implementation of a reversible function one can interchange controlled-V with 
controlled-y+ gates without changing the function realized, provided all such gates are interchanged. In our 
search procedure, this is accounted for by never using a controlled-y+ gate as the first quantum gate during 
construction of an NCV circuit. In this context, the first quantum gate is the one that transforms a Boolean 
line to one that can take on quantum values. 

• Once an optimal implementation of a function is found, we have also found an optimal implementation for 
all functions that differ from this one only by their input-output labeling. Potentially, this accounts for at 
most 6 different functions. 

• Once G\G2 - G]t is a circuit for a reversible function /, G^'G^Jj ...Gj"' is a valid circuit for 0. It 
can be shown that for each metric considered in this paper, as well as in any weighted linear type metric 
if GiG2 - Gic is optimal for /, then G^t'G^Jj ...Gj"' will be an optimal implementation for /^'. From the 
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point of view of the search for optimal circuits, this means that once an optimal circuit for / is found, so is 
an optimal circuit for /^'. This observation would further help to cut down the search space, however, we 
have not yet implemented it since our program is fast enough at present. It takes approximately 1 minute to 
synthesize optimal 3-line NCV circuits in each metric on a single 750 MHz processor Sun Blade 1000. 

Due to the potentially differing costs of the basis gates, our procedure maintains several queues of functions, 
each corresponding to the cost associated with the circuits in it. During the search, new gates are assigned to the 
circuits with smallest cost not yet considered thereby yielding new circuits to be considered. Cheaper gates are 
applied first. However, due to varying gate costs, the first circuit found realizing a Boolean function may not be 
optimal. To see this, consider the example illustrated in Figure^ Our program finds the non-optimal circuit with 
NCV-155 cost 7 before the optimal implementation with NCV-155 cost 6. This is because the procedure generates 
a circuit with two NOT gates before a circuit with a single CNOT gate, and consequently finds the three gate 
circuit in advance of the cheaper two gate alternative. Note that for gate count these two circuits are generated in 
the opposite order 
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(a) (b) 
Figure 1: (a) Non-optimal, but first found and (b) optimal in NCV-155 cost metric circuits. 



4 Observations on the optimal 3-line NCV circuits 

Using the above search procedure we found several interesting properties of the optimal 3-line NCV circuits: 

• Optimal implementations found with cost metrics NCV-11 1 and NCV-155 are interchangeable, and optimal 
NCV-111 implementations have optimal costs in NCV-012 cost metric. Diagram in Figure |2] shows which 
optimal implementations can be substituted without loosing the property of optimality in the corresponding 
metric. Practically, this means that the set of optimal NCV-111 circuits contains circuits optimal in other 
(NCV-012 and NCV-155) metrics. In Section|5lwe make some observations with regard to the NCV-111 
metric optimal circuits. The same comparisons apply for optimal NCV-012 and optimal NCV-155 imple- 
mentations. Further, our experiments with different metrics suggest that the set of optimal NCV-1 1 1 circuits 
will contain optimal implementations in NCV-xyz cost metric as long as non-negative integer numbers x, 
y, and z which represent costs of the gates NOT, CNOT, and controlled-V (and assuming the cost of the 
controlled-y+ equals the cost of the controlled-V) satisfy inequality y < 2z. 



Optimal X Optimal 

-* f NCV-111 J ^ * ( NCV-155 

xcircuits^/ \circuits 

Figure 2: Interchangeability of the optimal implementations. 



Optimal 
NCV-012 
circuits 



• In every circuit cost metric that we have considered (NCV-1 1 1, NCV-155, and NCV-012), the total number 
of controUed-y and controlled-y+ gates in any single circuit is divisible by 3 and is never more than 9. We 
conjecture that the overall number of controlled-V and controlled-y+ gates in any NCV implementation of 
any reversible function is divisible by three. 
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• A 3-bit Toffoli gate and a 3-bit Toffoli gate with one negative control have the same cost in each of the 
metrics NCV-111, NCV-155, and NCV-012. Some optimal circuits are illustrated in Figure|3] An optimal 
implementation of a 3-bit Toffoli gate with two negative controls has 6 gates. 
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Figure 3: (a) optimal NCV realization of the Toffoli gate (b),(c) optimal NCV realizations of the Toffoli gate with 
a single negative control. 



This observation allows us to generalize the famous result by Barenco et al. 121 and improve |6| on the 
implementation costs of the generalized multiple controlled Toffoli gates by suggesting that large Toffoli 
gates with some, but not all, negative controls can be simulated with the same cost as the same size Toffoli 
gate with all positive controls. Illustrated in Figurel^is a Toffoli gate with 5 controls, and a circuit simulating 
it, analogous to the one from f2\. Since the cost of a size 3 Toffoli gate with one negative control equals to 
the cost of size 3 Toffoli gate with both positive controls, large Toffoli gate illustrated in Figurel^will have a 
cost equal to the cost of the same size Toffoli gate but having only positive controls. Application of the local 
optimization techniques like (7) may choose different scenarios for simplification of the circuits for large 
Toffoli gates with all positive and some but not all negative controls. However, results of such optimization 
are out of the scope of the present paper. 
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Figure 4: Construction of a large Toffoli gate with some but not all negative controls. 



• Observing all optimal implementations in the NCV- 11 1 metric we came to the conclusion that every Boolean 
function f{x\ .xq) of two variables can be computed by a quantum circuit with no more than 5 gates. If this 
function is from Boolean class (when /(0,0) = 0), only 4 quantum gates are required. 

5 Comparisons of the Sets of Optimal Implementations 

In this section, we compare the NCV-111 (NCV-012 and NCV-155) costs for the optimal synthesis of size 3 
reversible functions using Toffoli gates up to size 3 (NCT library, 1141 ') with the costs of optimal NCV circuits. 
We note that our implementation of the breadth first search for the size 3 reversible NCT circuits may differ from 
the original as discussed in 1141 . Due to the large number of optimal NCT circuits for some functions, the results 
shown below may vary slightly depending on the actual program implementation. Further, the comparison is done 
through substitution of every Toffoli gate in an optimal NCT circuit with a 5-gate NCV circuit (111 |, page 182), 
and thus assigning a cost of 5 in the NCV-1 1 1 metric. Gates NOT and CNOT are present in both libraries, NCT 
and NCV, and thus require no specific attention. 
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Opt. NCT 


Opt. NCV-111 


Cost 


Gcm 


NCV-111 


NCV-111 





1 


1 


1 


1 


12 


9 


9 


2 


102 


51 


51 


3 


625 


187 


187 


4 


2780 


392 


417 


5 


8921 


475 


714 


6 


17049 


259 


1373 


7 


10253 


335 


3176 


8 


577 


1300 


4470 


9 





3037 


4122 


10 





3394 


10008 


11 





793 


5036 


12 





929 


1236 


13 





4009 


8340 


14 





8318 


1180 


15 





4385 





16 





255 





17 





1297 





18 





4626 





19 





4804 





20 





475 





21 





106 





22 





503 





23 





357 





24 





4 





27 





17 





28 





2 





WA 


5.8655 


14.0548 


10.0319 



Table 1: Optimal NCT and optimal NCV-1 1 1 NCV circuits. 



5.1 Optimal NCV-111 vs. Optimal NCT Circuits 

The leftmost column of Table [prefers to the cost of an implementation of a 3-bit reversible function. The second 
and third columns refer to the number of optimal NCT circuits reported in 1141 ; the fourth column refers to the 
number of optimal NCV circuits in the NCV-111 metric found by our search procedure. The circuit costs used 
in the last three columns are NCT gate count, NCV-111 and NCV-111 correspondingly. Table also^reports the 
weighted average (WA) for the given cost metric. The third and fourth columns compare the NCV-1 1 1 costs of the 
optimal NCT and optimal NCV-1 1 1 implementations in cost metric NCV-1 1 1 for the size 3 reversible functions. 
We observed that the maximal ratio of the cost of one optimal implementation over the other is 3.375 = ^. That 
is, even for circuits with a small number of inputs/outputs an optimal Toffoli circuit transformed to NCV can be 
a factor of 3.375 off its optimal NCV realization. We also determined that on average, the optimal NCT circuit 
is 1.3902 times more expensive (for the NCV-111 metric) than the corresponding optimal NCV circuit found by 
our search procedure. An important question of whether the NCT and NCV costs are related is addressed by the 
computing the correlation between the vectors of NCV-1 1 1 costs of the optimal NCT and optimal (in NCV-1 1 1 
metric) NCV circuits. The correlation equals 0.896. We conclude that the costs are reasonably correlated. Finally, 
we found it interesting to determine how many optimal NCT circuits have optimal NCV-1 1 1 cost. There are 1,610. 
The results are illustrated in the cost comparison chart in Figure|5] 
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Figure 5: Optimal NCT (X-coordinate) VS Optimal NCV-U 1 (Y-coordinate) circuits. 



5.2 Optimal NCV-012 vs. Optimal NCT Circuits 

We considered the set of optimal NCV-012 circuits as created by our program. Optimal NCV-11 1 implementations 
were not substituted for those optimal NCV-012 circuits whose NCV gate count is not optimal. 

A comparison of the NCV-012 costs of optimal NCV-012 NCV circuits and NCT circuits is made in Table |2] 
(this table is organized similarly to the previous table) and Figure |6l In this metric, the maximum ratio of NCT 
optimal circuit NCV-012 cost over NCV-012 optimal circuit NCV-012 cost equals 8 (= ^). The optimal NCT 
and optimal NCV-012 circuits for one such function are illustrated in Figure^. On average, however, this ratio 
is 1.2728. Correlation between the vectors of costs is 0.8999, which is almost identical to the correlation between 
optimal NCT and optimal NCV-111 circuits in the NCV-111 cost metric. The number of functions where the 
NCV-012 cost of an optimal NCT circuit equals the NCV-012 cost of an optimal NCV-012 circuit equals 1,774. 



6 Optimal implementations with restricted qubit-to-qubit interactions 

In the discussion above we assumed that a direct interaction between any two qubit can be established (e.g. a 
2-qubit gate may be built on any two qubits). However, due to the specifics of a particular physical realization, 
this may not always be the case. Some of the qubit-to-qubit interactions may only be available indirectly. On the 
other hand, in every «-bit quantum computation it must always be possible to construct a connected graph with 
vertices representing qubits and edges representing possibility of the direct interaction between qubits. In the case 
n = 3, there are only two non-isomorphic connected graphs. One is the complete graph, and the remaining is a 
star (all vertices are connected to one). In this section we report results for the optimal synthesis assuming direct 
interactions are allowed between qubits a and b, and b and c. We assume NCV-111 costing metric, however, the 
results can be calculated for other metrics as well. Assuming the implementation costs (numbers of gates) are 
0..23, the numbers of functions requiring this many gates are 1, 7, 29, 82, 181, 334, 374, 334, 337, 753, 1652, 
2654, 2482, 1674, 1350, 3236, 6304, 6028, 1508, 1302, 2566, 4314, 2804, and 14. There are no functions requiring 
more than 23 gates. 

We found it interesting to notice that in the case of non-restricted qubit interactions one of the cheapest non- 
linear (with respect to EXOR) reversible gates is the Peres gate 1121 defined by the transformation {a,b,c) 
{a,b®a,c®ab). It can be composed with 4 quantum NCV gates only. An analogy (smallest non-linear reversible 
gate with a similar transformation) of this gate in the case of restricted qubit interactions is the gate defined by the 
transformation {a,b,c) i-^ {b,a,c(Bab). Itrequires 6 quantumNCV gates as illustrated in FigurelSja). Toffoli gates 
TOF{a,b;c) and TOF{b,c;a) require 9 elementary quantum gates each and can be thought of as a SWAP gate 
(3 CNOTs) followed by the {a,b,c) i-^ {b,a,c ® ab) transformation. An optimal implementation of TOF{a,b;c) 
is illustrated in Figure|8jb). Toffoli gate TOF{a,c;b) is somewhat more expensive. It requires 13 elementary 
operations in its optimal implementation (see Figure |8jc)). 
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Opt. NCV-012 


Cost 


GC 


NCV-012 


NCV-012 


NCV-111 
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8 


8 


1 
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12 


48 


48 
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102 


183 


192 


45 
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625 


398 


408 


142 
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2780 


486 


480 


315 


5 
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201 


192 


585 
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17049 


16 


16 


1169 


7 


10253 





192 


2286 
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577 


47 


1056 


3414 


9 





352 


3168 


4790 


10 





1347 


4320 


6744 


11 





3130 


672 


6420 


12 





3340 





4328 


13 





561 





4360 


14 





3 


2880 


4032 


15 








11520 


1568 


16 





162 


4416 


112 


17 





1219 








18 





4435 








19 





8029 








20 





3872 








21 





128 


9856 





22 








896 





24 





341 








25 





1946 








26 





4482 








27 





3977 








28 





609 








29 





6 








32 





289 








33 





489 








34 





194 








35 





3 








40 





16 








41 





3 








WA 


5.8655 


19.1011 


14.9800 


10.5800 



Table 2: Optimal size 3 reversible circuit NCV-012 costs in NCT and NCV bases. 
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Figure 6: Optimal NCT (X-coordinate) VS Optimal NCV-012 (Y-coordinate) circuits. 
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Figure 7: (a) optimal NCT and (b) optimal NCV-012 circuits for the function [7,6,4,5,2,3, 1,0]. 



7 Conclusion 

The results in Section|5]lead us to the following conclusion. Minimization of Toffoli gate count as a criterion for 
a reversible synthesis method is not optimal and even for small parameters may result in a seemingly small circuit 
which may be as far off a technologically favorable implementation as a factor of 8. It is natural to expect that for 
circuits with more lines the difference will grow. We suggest that the commonly used gate count metric should be 
replaced with a metric that accounts for the different costs of large building block (i.e., Toffoli gates), such as the 
weighted gate cost presented here. Using such a metric would lead to the technologically favorable circuit (b) in 
Figure 

Finally, we observe that for some small parameters minimizing the cost in any of the technology oriented 
metrics that we considered in this paper (NCV-1 1 1, NCV-155 or NCV-012) should result in a reasonably small if 



(c) 







(a) 






f 9 t 

i 


— — T — ■ 








y CD * _T_ CD y 








■a — 01 





Figure 8: Optimal circuits for (a) {a,b,c) i-^ {b,a,c(Bab), (b) TOF{a,b;c), and (c) TOF{a,c;b). 
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not optimal solution if another metric is applied to the cascade of gates. Thus, at the present time it is sufficient to 
use any of the suggested technology oriented cost metrics. NCV-111 and NCV-155 are preferred over the NCV- 
012 metric, since all such reahzations (assuming reversible functions with no more than 3 variables) have optimal 
NCV-012 cost. 
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